Hansen et al. Journal of Biological Engineering 2012, 6:1 1 
http://www.jbioleng.0rg/content/6/l/l 1 



JOURNAL OF BIOLOGICAL 
ENGINEERING 



METHODOLOGY 



Open Access 



Development of a modified Respiration Activity 
Monitoring System for accurate and highly 
resolved measurement of respiration activity in 
shake flask fermentations 

Sven Hansen, loanna Hariskos, Bettina Luchterhand and Jochen Buchs* 



Abstract 

Background: The Respiration Activity Monitoring System (RAMOS) is an established device to measure on-line the 
oxygen transfer rate (OTR), thereby, yielding relevant information about metabolic activities of microorganisms and 
cells during shake flask fermentations. For very fast-growing microbes, however, the RAMOS technique provides too 
few data points for the OTR. Thus, this current study presents a new model based evaluation method for 
generating much more data points to enhance the information content and the precision of OTR measurements. 

Results: In cultivations with E.coli BL21 pRSET eYFP-IL6, short diauxic and even triauxic metabolic activities were 
detected with much more detail compared to the conventional evaluation method. The decline of the OTR during 
the stop phases during oxygen limitations, which occur when the inlet and outlet valves of the RAMOS flask were 
closed for calibrating the oxygen sensor, were also detected. These declines reflected a reduced oxygen transfer 
due to the stop phases. In contrast to the conventional calculation method the new method was almost 
independent from the number of stop phases chosen in the experiments. 

Conclusions: This new model based evaluation method unveils new peaks of metabolic activity which otherwise 
would not have been resolved by the conventional RAMOS evaluation method. The new method yields 
substantially more OTR data points, thereby, enhancing the information content and the precision of the OTR 
measurements. Furthermore, oxygen limitations can be detected by a decrease of the OTR during the stop phases. 

Keywords: Oxygen transfer rate, RAMOS, Shake flask, Bioprocess monitoring, Diauxic growth, Oxygen limitation 



Background 

Shake flasks are widely used in fermentations for biotech- 
nological research and industrial process development [1,2]. 
For gaining a better understanding and control of shake 
flask cultivations, various methods have been recently 
developed for online monitoring of process parameters. 

Monitoring of pH-values in shake flasks has been rea- 
lized both with standard autoclavable pH-probes that are 
immersed in the bulk liquid [3] and with fluorescent 
optodes fixed at the flask wall that allow optical measure- 
ment [4-6]. Moreover, dissolved oxygen tension (DOT) 
can in principle be measured in shake flasks by using 
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either Clark- type electrodes [7-10] or optical sensors based 
on dynamic quenching of luminescence [11-16]. These 
non-invasive measurement methods have proven to be 
more reliable, since they do not alter the hydrodynamics 
of the culture system due to baffling effects [17,18]. 

As shown in numerous studies, almost all metabolic 
activities of aerobic microorganisms depend on the oxygen 
consumption of the culture [8,19-21], Therefore, online 
measuring techniques are also useful for determining the 
gas transfer rates in shake flasks. The company BlueSens 
GmbH (Herten, Germany) has developed a system that 
measures the gas transfer rates in the headspace of shake 
flasks down to nominal flask volumes of 500 mL [22]. 

As an alternative method, the Respiration Activity 
Monitoring System (RAMOS) presented by Anderlei 
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et al [23,24], represents another non-invasive measure- 
ment technique which allows the online determination of 
oxygen transfer rate (OTR), carbon dioxide transfer rate 
(CTR) and respiratory quotient (RQ) in shake flasks down 
to 100 ml Since its introduction, RAMOS has been used 
for multiple applications: e.g., for determining oxygen lim- 
itations in shake flask cultivations [18,25-30], screening of 
microorganisms [31-33], optimizing media [34-36], inves- 
tigating secondary substrate limitations [37,38] and stress 
phenomena [39], process development and optimization 
[2,40-44] and for monitoring precultures for fermentations 
in stirred tank reactors [45,46]. The setup of a RAMOS 
device is illustrated in Figure 1. For determining OTR with 
RAMOS, a certain measuring cycle is regularly repeated 
and consists of a rinsing phase and a stop phase [23,24]. 
In the rinsing phase air is injected into the RAMOS flask 
to ensure a headspace gas concentration that equals that 
in a normal Erlenmeyer flask [24]. During the stop phase 
the air flow into the system is briefly interrupted, and then 
the OTR is calculated from the corresponding drop (linear 
slope) in the oxygen partial pressure as detected by an 
oxygen sensor. Since these measuring cycles are typically 
repeated every 30 min, a data density of merely two OTR 
measurement points per hour is usually achieved. 

The objective of this study is to present a new model 
based evaluation method for generating substantially more 
data points to enhance the information content and the 
precision of OTR measurements. First, the conventional 
way of generating OTR data is recalled. Then, the new 
evaluation method based on a complete oxygen headspace 
balance is developed. Microbial cultivations have been 
devised in order to evaluate the new method. 

Results and discussion 

Conventional evaluation method of the RAMOS device 

In the conventional approach as proposed by Anderlei 
et al. [23,24], the OTR in RAMOS is calculated during 



the stop phase based on equation A.l. The oxygen sen- 
sor does not directly provide the oxygen partial pressure 
P02 but a voltage signal U Q 2 that depends on it. There- 



fore, the expression 



Po 2 , 



Uo 2 . 



is introduced to convert this 



voltage signal into the oxygen partial pressure during the 
stop phase as described by the right part of equation 
A.l. Here, V G and V L are the volumes of the gas head- 
space and the liquid in the flask, respectively, and T is 
the temperature. 
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To compensate for the signal drift, the sensors are 
calibrated by calculating the partial pressure of oxygen 
at the end of the rinsing phase po 2 ,reai by applying the 
steady state gas composition of the headspace volume at 
0 h. Figure 2 gives an example of the sensor signal dur- 
ing a stop phase. This leads to equation A.2 for calculat- 
ing the OTR: 
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^top 
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Here, po2,in is the oxygen partial pressure of the inlet 
flow, and V in is the volumetric flow into the flask during 
the rinsing phase, p am b is the ambient pressure and V m 
is the molar volume. 

A calibration factor (K) can then be calculated by 
using equation (A.3). 
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Figure 1 Setup of a RAMOS device as introduced by Anderlei et al. [23,24]. 
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Figure 2 Example of a calibration during a RAMOS measurement, with measured data from the oxygen sensor {circle), smoothing 
spline fit of the rinsing phase according to equation (B.3) {green line), optimized polynomial of the stop phase according to equation 
(B.5) used for solving equation (B.4) {blue line) and delayed signal of the stop phase according to the delay model described by 
equation (B.6) {red line). 
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To examine the conventional evaluation method pre- 
sented by Anderlei et al [23,24] Exoli BL21 pRSET 
eYFP-IL6 was cultivated in Wilms-MOPS-medium as 



described in "Materials and methods". Figure 3A illus- 
trates the OTRs calculated using equation (A.2). Here, 
the OTR curve consists of one measuring point every 
30 min, because the OTR is only evaluated during the 
stop phases. This curve shows a typical diauxic growth 
with the first peak leveling off and forming a horizontal 
plateau between 9.5 h and 10.5 h after inoculation. At 




Time [h] 

Figure 3 Oxygen transfer rates (OTR) and calibration factors (K) of cultures of E. coli BL21 pRSET eYFP-IL6 in mineral medium 
measured with the RAMOS device by applying both the method of Anderlei et al. [23,24] and the newly developed evaluation 
method. A: OTR calculated by applying the evaluation method of Anderlei et al. [23,24] using equation A.2 {block triangle) and the corresponding 
calibration factor K calculated with equation A3 {red triangle); B: OTR calculated by applying the new evaluation method using equation (B.1) in 
the rinsing phase {open square), equation (B.4) in the stop phase {black square) and K, calculated as illustrated in Figure 2 and equation (B.7) 
{red square) with its smoothing spline fit {red line); experimental conditions: Wilms-MOPS medium with 20 g/L glucose; temperature T= 30°C; 
filling volume V L = 1 0 ml_; shaking diameter d 0 = 50 mm; shaking frequency n = 350 rpm; initial optical density OD 600 = 0.2; initial pH = 7.5; 
rinsing phase t rinse = 25 min; stop phase t stop = 5 min; sensor lag time t = 0.01 3 h. 
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the end of this phase glucose is fully consumed under 
oxygen limitation [23,24]. At around 12.5 h after inocu- 
lation, a second peak appears indicating the consump- 
tion of acetate formed during the first phase of glucose 
consumption. This typical growth behavior has already 
been confirmed in other studies [6,47,48]. Therefore, the 
conventional evaluation method of RAMOS is able to 
measure the respiration activity of shake flask cultures and 
to detect phenomena as oxygen limitations and diauxic 
growth behavior, as long as these effects are lasting longer 
than the typical measuring cycle of 30 min. This condition 
is valid for many biological systems. However, a higher 
data density would be desirable when a highly resolved 
measurement is needed, e.g. for tracking metabolic activ- 
ities that take place within much shorter periods. 

Also illustrated in Figure 3A is the calibration factor K as 
calculated by equation (A.3). This curve also shows two 
peaks which both occur right after the OTR peaks. These 
two peaks are totally unexpected, because the calibration 
factor K is assumed to be steady and independent from the 
OTR. Consequently, the assumption that the gas headspace 
volume at the end of the rinsing phase is in a steady state 
results in an inadequate calibration. Therefore, a calculation 
method should be developed that increases the data density 
of the OTR and includes a more accurate calibration. 

New model based evaluation method 

For increasing the data density of the OTR measurement 
the data from the rinsing phase should also be utilized 
for the evaluation. Therefore, the complete oxygen bal- 
ance of the headspace of a RAMOS flask is derived and 
resolved for the OTR as follows (equation (B.l)): 



OTR ri 



V in 'Po2M - V out 'K'Uo 2 , 
V L RT 

dUo 2 .rinse Vg 



K 



dt V L -R-T 



(B.l) 



V in and V out describe the volumetric flows into and 
out of the flask, respectively, po2,in the oxygen partial 
pressure of the inlet flow, K is the calibration factor of 
the oxygen sensor, U Q 2 is the oxygen sensor signal, V G 
and V L are the volumes of the gas headspace and the 
liquid medium in the flask, respectively and T is the 
temperature. 

The outlet flow V out is calculated from the slope of the 
pressure sensor signal of the headspace volume in the 
adjacent stop phases according to equation (B.2). This 
change in pressure during the stop phases occurs when 
the respiratory quotient (RQ) differs from 1. During the 
rinsing phases, when the valves are opened, this results in 
a change of the outlet flow. Equation (B.2) also takes into 
account the water vapor partial pressure that affects the 
outlet flow. 
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Here, p am b describes the ambient pressure and Ph20 is 
the water vapor partial pressure in the gas headspace 
volume. 

For determining the slope of the oxygen sensor signal 
dU °^ inse , the sensor data of each rinsing phase are 
approximated with a smoothing spline function f that 
minimizes equation (B.3) according to the algorithm of 
de Boor [49]. This algorithm is implemented into the 
Curve Fitting Toolbox of MATLAB R2010b (The Math 
Works, MA, USA). The left part of equation (B.3) con- 
sists of a least square method and is responsible for a 
close approximation of the sensor signal. In contrast, the 
right part of the equation represents the second deriva- 
tive of the fit and, therefore, is responsible for the 
smoothness required for calculating reliable derivations. 
For adjusting a good trade-off between these two terms 
the regularization parameter p is used. A regularization 
parameter of p = 0.999 is suitable to describe the sensor 
signal of the rinsing phases and is, therefore, used to 
perform the fitting 
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p is the regularization parameter for the smoothing 
spline fitting and f describes the smoothing spline function. 

For calculating the OTR in the stop phase, equation 
(B.l) simplifies to equation (B.4) since no air flow is 

applied to the flask. Therefore, V in and V out are zero 
and the left part of equation (B.l) vanishes. 
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The oxygen partial pressure in the stop phase is approxi- 
mated with a second-order polynomial as described by 
equation (B.5). The reason for using a second-order poly- 
nomial equation is because it also considers a change of 
OTR during the stop phase and uses just one parameter 
more than a linear slope. As the stop phases are relatively 
short a more complicated equation resulting in a higher 
computational effort does not seem to be reasonable. 
Since the changes of the OTR are not restricted to the rin- 
sing phase this method can be considered to be more 
accurate than using a linear slope. The first derivative of 
this polynomial is inserted into equation (B.4). 
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Due to the inertia of the oxygen sensor, there is a lag 
time for the sensor signal which does not allow a direct fit 
of equation (B.5). This lag time is compensated for by 
incorporating the first order delay model equation (B.6) 
into the fitting procedure. Thereby, the calculation of the 
coefficients a, b and c of the polynomial (B.5) is based on a 
least-square optimization procedure of these coefficients. 
Equation (B.6) is used to constantly calculate the corre- 
sponding sensor response U 0 2,deiay on equation (B.5) by 
varying the estimates of a, b and c. As soon as the devia- 
tion between the calculated sensor response U 0 2,deiay and 
the measured signal is small enough, the optimization 
stops and the course of the oxygen partial pressure can be 
calculated by using the corresponding estimates of the 
coefficients for the polynomial (B.5). Figure 2 also shows 
the optimized polynomial of the stop phase according to 
equation (B.5) and the delayed signal of the stop phase 
according to the delay model described by equation (B.6). 
For the oxygen sensors used in this work, a time constant 
t of 0.013 h is used (see Materials and methods). This value 
best describes the dynamic behavior of the sensor. 

dUo 2j delay Uo 2 ,stop ~ Uo 2 , delay (v> r\ 

~^— = r (B - 6) 

Sensor calibration 

For solving equations (B.l) and (B.4), the calibration factor 
K is required. A new calibration method has been devel- 
oped which also considers the dynamic behavior of the gas 
headspace volume instead of assuming a steady state. The 
end of the rinsing phase as well as the start of the subse- 
quent stop phase of a biological experiment is utilized in 
the new sensor calibration method. This point occurs at 
0 h in Figure 2. 

Figure 2 also shows a section of the sensor signal curve 
during a RAMOS measurement and the mathematical fits 
including both the rinsing phase and the stop phase. 
Within the rinsing phase which last up to the beginning of 
the stop phase at 0 min, the sensor signal is described by a 
smoothing spline function f according to equation (B.3). 
During the stop phase, from 0 min to 5 min, the sensor 
signal strongly decreases. However, due to the lag time at 
the beginning of this stop phase, the initial slope of the 
signal curve merely changes slowly. Figure 2, thus, shows 
both the polynomial described by equation (B.5) and the 
delayed signal described by equation (B.6) after the opti- 
mization procedure. 

To calculate the calibration factor K, equations (B.l) 
and (B.4) are used. Whereas equation (B.l) describes the 
OTR during the rinsing phase and equation (B.4) the 
OTR during the stop phase, both equations can be 
equalized at the intersection of the curves of the rinsing 
phase and the stop phase as illustrated in Figure 2. This 



is justified, because a sudden change in the OTR or K 
due to a change in the gas flow can be excluded. By 
resolving equations (B.l) and (B.4) under the assumption 
of their equality (OTR rinse = OTR stop ), equation (B.7) is 
obtained to calculate K. 

K = - V in' P Q2,in 

/ dllQ2^stop, inter sec dUo2,rinse, inter sec \ ~ir ~ir j r 

\ dt dt ) ' Vg ~ V out' U 0 2 , inter sec 

(B.7) 

Using the same experiment, Figure 3B depicts the 
OTR curve resulting from the new calculation method. 
Here, an OTR is calculated every 5 min. The OTR data 
from the stop phase are in good agreement with the 
OTR data calculated using the method of Anderlei et al. 
[23,24]. The OTR data from the rinsing phase mostly 
concur with that of the stop phase. Only during the 
phase of oxygen limitation, indicated by the horizontal 
plateau of the OTR data at ca. 9.5 - 10.5 h [23,24], the 
OTR in the stop phase is somewhat lower than that in 
the rinsing phase. 

To understand this phenomenon, Figure 4 is considered. 
The oxygen partial pressure of the gas headspace volume 
and the dissolved oxygen concentration in the medium of 
a simulation example are shown. The simulation is based 
on the model of equations (C.l) to (C.7) as shown in the 
appendix. This simulation example describes an oxygen 
limitation between 19 h and 21 h. The dissolved oxygen 
concentration c Q2 was calculated by using the oxygen 
transfer equation (B.8). A constant mass transfer coeffi- 
cient k L a was assumed. 

OTR = k L a.(c 02 * -c 02 ) (B.8) 

Up to the beginning of the oxygen limitation at 19 h 
both the oxygen partial pressure of the headspace and the 
dissolved oxygen generally decrease due to increasing 
respiration activity. However, during the rinsing phases of 
this oxygen limitation between 19 h and 21 h the oxygen 
partial pressure of the headspace shows a horizontal pla- 
teau of ca. 0.18 bar whereas the dissolved oxygen has very 
low values which are close to the K Q 2 value of oxygen. 
Both the oxygen partial pressure of the headspace and the 
dissolved oxygen generally increase after the oxygen lim- 
itation when the respiration activity decreases. 

Up to the oxygen limitation at 19 h, both, the oxygen 
partial pressure of the gas headspace volume and the dis- 
solved oxygen strongly decrease during the stop phases as 
can be seen by the spikes in the signal. As both signals 
decrease the driving force of the oxygen transfer into the 
medium is not affected. However, during oxygen limita- 
tion between 19 h and 21 h the dissolved oxygen in the 
medium approaches very low values in the range of the 
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Time [h] 

Figure 4 Oxygen partial pressure of the gas headspace volume and dissolved oxygen during oxygen limitation. Example of a RAMOS 

fermentation with oxygen partial pressure of the gas headspace volume {block line) and dissolved oxygen in the medium (red line) calculated 

with the model of equations (C.1) to (C.7) in the appendix, 
k J 



K M value for oxygen and cannot decrease further. Never- 
theless, the oxygen partial pressure of the headspace 
decreases, resulting in a lower driving force for the oxygen 
transfer into the medium and, consequently, in a slightly 
lower OTR during the stop phases, as depicted for 9.5 h- 
10.5 h in Figure 3B. 

In Figure 3B also the calibration factor K as calculated 
with equation (B.7) is shown. In contrast to the steady 
state calibration using equation (A.3) as shown in 
Figure 3A, the calibration factor K calculated with equa- 
tion (B.7) has a steady course as expected and does not 
show any peaks correlated to the OTR. However, the 
data for K show some fluctuations which are due to the 
lower accuracy of data fitting at higher sensor dynamics. 
Therefore, the data for K are smoothed by using a 
smoothing spline with a regularization parameter of 
p = 0.1. This value has been found to be a good trade-off 
between smoothing of K and revealing the effect of 
sensor drift. 

Case studies 

Recombinant E. coli fermentation in mineral medium 

To verify the new evaluation method, E.coli BL21 pRSET 
eYFP-IL6 was cultivated using a RAMOS device with 
Wilms-MOPS -medium containing glucose and sorbitol 
as carbon sources. Four shake flasks were examined, 
each having a different initial sorbitol concentration. 
Figure 5 shows the results. According to Monod [50], 
this cultivation leads to a diauxic growth where first glu- 
cose is consumed followed by a consumption of sorbitol. 
In all four cases, the OTR curves show an exponential 
growth phase during the 8 h after inoculation. After- 
wards, the OTR curves form a horizontal plateau for up 



to 9.75 h after inoculation which indicates an oxygen 
limitation during this period of time [23,24]. This oxygen 
limitation becomes more apparent when considering 
that the OTR of the stop phases are slightly lower during 
this period compared to the OTR of the rinsing phases. 
This difference is attributed to a reduced driving force of 
the oxygen transfer during the stop phases at oxygen 
limitation as described above. The oxygen concentration 
in the medium approaches very low values in the range 
of the K02 value for oxygen, whereas the oxygen partial 
pressure of the headspace is still decreasing. Up to 
9.75 h after inoculation when the oxygen limitation ends 
and the OTR decreases, glucose is completely con- 
sumed. Before glucose is depleted, the sorbitol concen- 
trations remain constant and acetate is formed as an 
overflow metabolite of E. coli fermentation [6,47,48]. 
Furthermore, the optical density increases indicating 
biomass formation. 

After 9.75 h, when the glucose is fully exhausted, sorbi- 
tol is consumed as can be seen by a decline in the sorbi- 
tol concentration. This sorbitol consumption is also 
observed in the OTR curve via a peak beginning at 10 h. 
The higher the initial concentration of sorbitol, the more 
pronounced the peak of the OTR curve is (from 
Figure 5 A to 5D). Whereas the second OTR peak reaches 
a maximum of 0.039 mol/L/h at an initial sorbitol con- 
centration of 0.5 g/L (Figure 5A), it attains a maximum 
of 0.048 mol/L/h at an initial sorbitol concentration of 
1 g/L (Figure 5B) and a limiting value of 0.063 mol/L/h 
at an initial sorbitol concentration of 1.5 g/L (Figure 5C). 
At an initial sorbitol concentration of 2 g/L (Figure 5D), 
the OTR peak also reaches limiting values at about 
0.062 mol/L/h, but is much wider. Thus, due to the 
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Figure 5 Oxygen transfer rates (OTR) of E. coli BL21 pRSET eYFP-IL6 in mineral medium with different initial sorbitol concentrations. 

A: 0.5 g/L sorbitol; B: 1 g/L sorbitol; C: 1,5 g/L sorbitol; D: 2 g/L sorbitol; RAMOS data from the rinsing phase calculated with equation (B.1) 
{open square) and from the stop phase calculated with equation (B.4) {black square); glucose concentration {blue triangle); sorbitol concentration 
{orange circle); acetate concentration {green diamond); optical density (OD 600 ) at 600 nm {red open circle); experimental conditions: Wilms-MOPS 
medium with 20 g/L glucose and different sorbitol concentrations; temperature T = 30°C; filling volume V L = 1 0 ml_; shaking diameter d 0 = 50 mm; 
shaking frequency n = 350 rpm; initial optical density OD 600 = 0.2; initial pH = 7.5; rinsing phase t rinse = 25 min; stop phase t stop = 5 min; sensor lag 
time t = 0.01 3 h. 



higher data density, the new evaluation method unveils 
the real size of these peaks. 

When sorbitol is depleted, the OTR curve again 
decreases before rising up to form a third peak which 
infers that acetate is now being consumed. This third peak 
has the same size for all four initial sorbitol concentrations, 
because the acetate is formed during the phase of glucose 
consumption and not during the sorbitol consumption. 
When acetate is consumed biomass is not produced any- 
more. Consequently, in the historical work of Monod [50], 
who measured biomass, only diauxic growth was observed 
and not a triauxic metabolic activity, as observed with the 
newly proposed calculation method of RAMOS data. 



To show the strengths of the new evaluation method, 
three fermentation examples were considered where the 
OTR data of both methods are shown. Figure 6A illustrates 
an enlarged section of the OTR curve of Figure 5A. This 
OTR curve indicates an oxygen limitation between 8.5 h 
and 9.5 h as the OTR data points of the rinsing phases 
form a horizontal plateau during this time [23,24]. Subse- 
quently, the OTR drops sharply, indicating a reduced 
respiration activity until 11 h after inoculation. The OTR 
curve calculated with the new method generates a mea- 
surement point every 5 min, hence, offers a very high data 
density. Therefore, a second peak between 10 h and 10.5 h 
could be detected which is due to the consumption of 
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Figure 6 Comparison of the method of Anderlei et al. [23,24] using equation (A.2) {blue triangle) and the new evaluation method using 
equation (B.I) in the rinsing phase (open square) and equation (B.4) in the stop phase {black square). A: Enlarged section of the OTR 
curve of Figure 5A; B: E. coll BL21 pRSET eYFP-IL6 in TB medium; temperature T=30°C; filling volume V L = 10 ml_; shaking diameter d 0 = 50 mm; 
shaking frequency n = 300 rpm; inoculated with 0.2 ml cryo culture; rinsing phase t rinse = 55 min; stop phase t st0 p = 5 min; high flow phase 
thigh = 0.9 min; C: Gluconobacter oxydans 621 H wildtype in complex medium with 40 g/L mannitol and 0.5 g/L glycerol; temperature T= 30°C; 
filling volume V L = 1 0 ml_; shaking diameter d 0 = 50 mm; shaking frequency n = 350 rpm; initial optical density OD 600 = 0.1 ; rinsing phase 
t rinse = 25 min; stop phase t st0 p = 5 min; a sensor lag time of t = 0.01 3 h is used for all experiments. 



sorbitol as could be shown in Figure 5A. This shows that 
the newly proposed method is very useful in displaying 
short-term effects of the OTR. 

The OTRs calculated with the new evaluation method 
in the stop phase are slightly higher than those calcu- 
lated with the method of Anderlei et al. [23,24]. This 
deviation is caused by the consideration of the sensor 
lag time in the new calculation method (equation (B.6)). 

Using the new method, it is clearly observed that during 
the oxygen limitation between 8.5 h and 9.5 h the data 
points in the stop phase (filled squares) are not in line with 
the data points from the rinsing phase (open squares), in 
contrast to the rest of the fermentation. As discussed pre- 
viously, this difference is due to a reduced driving force of 
the oxygen transfer during the stop phases of an oxygen 
limitation. The oxygen concentration in the medium 



approaches very low values, whereas the partial pressure 
of oxygen in the gas headspace volume still decreases. This 
effect does not occur in usual Erlenmeyer flasks, where 
the diffusive mass transfer through the cotton plug is not 
interrupted. Consequently, the OTR of the stop phase, 
gives a slightly different OTR than appearing in a typical 
Erlenmeyer flask. 

Recombinant E. coli fermentation in complex medium 

Figure 6B illustrates the OTR curves of a RAMOS culti- 
vation of Exoli BL21 pRSET eYFP-IL6 in TB-medium. 
Instead of the typically chosen value of 25 min, the rin- 
sing phase was set to a duration of 55 min. However, the 
data density of the OTR curve calculated with the new 
evaluation method is not affected by the reduced num- 
ber of stop phases because the data of the rinsing phase 
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are also considered. After the lag phase at 7 hours the 
OTR shows an exponential growth up to 10.5 h leading 
directly into a horizontal plateau which can be considered 
as an oxygen limitation [23,24] . Afterwards, even two very 
distinct peaks can be seen between 13 h and 16 h before 
the OTR drops down to a level of 0.005 mol/L/h at ca. 
17 h. This clearly shows that the new method is also inde- 
pendent from the number of stop phases selected in the 
experiments. This characteristic is advantageous in appli- 
cations where short-term effects are expected but a high 
number of stop phases is not desired, e.g. for preventing 
further oxygen limitations due to frequent interruption of 
the air flow. 

Gluconobacter oxydans fermentation in complex medium 

Figure 6C illustrates the OTR curves of an RAMOS cultiva- 
tion of Gluconobacter oxydans 621 H wild type in complex 
medium with 40 g/L mannitol. A rinsing phase of 25 min 
was selected. The OTR curve calculated with the method 
of Anderlei et al. [23,24] indicates an exponential growth of 
the culture up to 7 h of cultivation when an OTR of 
0.035 mol/L/h is achieved. The next measuring point at 
7.5 h also indicates an OTR of 0.035 mol/L/h before the 
OTR decreases again. Consequently, the time between 7 h 
and 7.5 h could be interpreted as a short period of oxygen 
limitation. Even though the OTR curve calculated with the 
new method basically shows the same course, the time span 
between 7 h and 7.5 h shows considerably more data 
points, and more importantly, the OTR curve with the new 
calculation method does not show an oxygen limitation 
between the two stop phases. 

Conclusions 

The newly proposed evaluation method yields substantially 
more OTR data points than the conventional method by 
Anderlei et al. [23,24]. This new evaluation method unveils 
additional peaks of metabolic activity which otherwise 
would remain undetected by the former method. Conse- 
quently, this new technique is a sophisticated means to 
generate more detailed information about metabolic activ- 
ities of any kind of microorganisms and cells during shake 
flask cultivations. Additionally, possible oxygen limitations 
can be detected by a decrease of the OTR during the stop 
phases of the RAMOS measurement. 

Materials and methods 

Organisms 

Exoli BL21 pRSET eYFP-IL6 was maintained at -80°C in 
Lysogeny broth (LB) medium with 100 ug/mL ampicillin. 
Stock solutions contained 200 g/L glycerol. Gluconobacter 
oxydans 621 H wild type was maintained in its cultivation 
medium (see below) including 80 g/L mannitol at -80°C. 
Stock solutions contained 150 g/L glycerol. 



Media 

LB medium for maintaining Exoli consists of: 5 g/L yeast 
extract (powdered, Roth, Karlsruhe, Germany), 10 g/L tryp- 
tone (pancreatic digest of casein, Roth, Karlsruhe, Germany) 
and 10 g/L NaCl. Terrific broth (TB) medium was used for 
cultivating both Exoli precultures and main cultures. The 
medium consists of: 5 g/L glycerol, 12 g/L tryptone (pan- 
creatic digest of casein, Roth, Karlsruhe, Germany), 24 g/L 
yeast extract (powdered, Roth, Karlsruhe, Germany), 
12.54 g/L K 2 HP0 4 and 2.31 g/L KH 2 P0 4 . Additionally, 
0.1 g/L ampicillin was added. 

Main cultures of Exoli in mineral medium were culti- 
vated in modified Wilms & Reuss synthetic medium (hen- 
ceforth referred to as Wilms-MOPS medium) [6,51]. This 
medium consists of: 20 g/L glucose, 5 g/L (NH 4 ) 2 S0 4 , 
0.5 g/L NH 4 C1, 3 g/L K 2 HP0 4 , 2 g/L Na 2 S0 4 , 0.5 g/L 
MgS0 4 7H 2 0, 41.85 g/L (0.2 M) 3-(N-morpholino)-propa- 
nesulfonic acid (MOPS), 0.01 g/L thiamine hydrochloride, 
1 mL/L trace element solution (0.54 g/L ZnS0 4 7H 2 0, 
0.48 g/L CuS0 4 -5H 2 0, 0.3 g/L MnS0 4 -H 2 0, 0.54 g/L 
CoCl r 6H 2 0, 41.76 g/L FeCl 3 -6H 2 0, 1.98 g/L CaCl r 2H 2 0, 
33.39 g/L Na 2 EDTA (Titriplex III)) and concentrations of 
sorbitol ranging from 0-2 g/L. Additionally, 0.1 g/L ampi- 
cillin was added. 

Gluconobacter oxydans was cultivated in a complex 
medium consisting of: 5 g/L yeast extract (powdered, 
Roth, Karlsruhe, Germany), 1 g/L K 2 HP0 4 , 1 g/L (NH 4 ) 
2 S0 4 , 0.5 g/L glycerol, 2.5 g/L MgS0 4 7H 2 0 and 40 g/L 
mannitol. The pH- value was adjusted to 6 with 1 M KOH. 

Cultivation conditions 

E. coli precultures and main cultures were cultivated in 
250 ml shake flasks at a temperature of T = 30°C with a 
shaking diameter of d 0 = 50 mm and a shaking frequency 
of n = 350 rpm (shaking machine LS-W, Kuehner AG, 
Birsfelden, Switzerland). For E. coli precultures, 20 ml of 
TB medium was inoculated with 200 ul of stock solution. 
Main cultures of Exoli in 10 mL Wilms-MOPS medium 
were inoculated with preculture broth from the exponential 
growth phase, resulting in an optical density (OD) of 0.2. 

The main culture of Exoli in TB medium was cul- 
tivated in 250 ml shake flasks with a filling volume of 
V L = 10 mL, a shaking diameter of d 0 = 50 mm, a shak- 
ing frequency of n = 300 rpm, and at a temperature of 
T = 30°C. This culture was inoculated with 200 ul of 
stock solution. 

The preculture and the main culture of Gluconobacter 
oxydans were cultivated in 250 ml shake flasks with a 
filling volume of V L = 10 mL, a shaking diameter of 
d 0 = 50 mm, a shaking frequency of n = 350 rpm, and at 
a temperature of T = 30°C. The preculture was inocu- 
lated with 500 uL of stock solution. The main culture 
was inoculated with preculture broth from the exponen- 
tial growth phase resulting in an OD of 0.1. 
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Respiration Activity Monitoring System (RAMOS) 

The respiration activity was measured in modified 250 mL 
Erlenmeyer flasks using a self-made RAMOS device as 
introduced by Anderlei et al [23,24]. The setup of the 
RAMOS device is shown in Figure 1. The 8 RAMOS flasks 
in parallel were supplied with air using a mass flow con- 
troller (type 5850 TR, Brooks, Hatfield, PA, USA). The air 
flow through the flasks is controlled by using valves at the 
inlet and the outlet of each RAMOS flask. For the rinsing 
phase the RAMOS flasks were flushed for 25 min with air 
at a flow rate of (low flow), which corresponds to an aera- 
tion rate of 1 wm at a filling volume of V L = 10 mL. The 
rinsing phase was followed by a 5 min stop phase (t stop ) 
with no air flow through the flasks. To compensate for the 
drop in oxygen partial pressure during the stop phase, a 

higher air flow of V in high = 3.6L N /h (high flow) was 
applied during the first 0.9 min of the rinsing phase. This 
measuring cycle is regularly repeated. The oxygen partial 
pressure of the headspace in a RAMOS flask was measured 
with a MAX250 oxygen sensor from Maxtec (Salt Lake 
City, Utah, USA). The difference between the headspace 
pressure and ambient pressure was detected with a 26PCA 
pressure sensor from Honeywell (Morristown, NJ, USA). 
The RAMOS flasks ensure the same hydrodynamic condi- 
tions and headspace gas concentrations as are found in 
regular Erlenmeyer flasks with cotton plugs [24]. Commer- 
cial versions of the RAMOS device are available from 
Kuehner AG, Birsfelden, Switzerland and Hitech Zang, 
Herzogenrath, Germany. 

The sensor lag time x has been determined in step 
change experiments. Beginning at a steady state of air in 
the gas headspace volume, oxygen depleted air was 
flushed into a RAMOS flask under well-defined condi- 
tions (y in = 0.54 mL/min, p 02 ,in = 0.197 bar, V L = 10 mL, 
n = 300 rpm). By calculating the real partial pressure of 
oxygen in the gas headspace of the flask and then opti- 
mizing the sensor lag time x, the model described by 
equation (B. 6) was fitted to the oxygen partial pressure 
signal measured by the sensor. This leads to a sensor lag 
time of x = 0.013 h (data not shown). 

Parallel shake flask cultivations 

Samples were taken from E. coli cultivations in Wilms- 
MOPS medium in Erlenmeyer flasks in parallel to the 
RAMOS experiments and cultivated under the same condi- 
tions as used in the RAMOS experiments. The OD of the 
samples at 600 nm was determined with a Thermo Scienti- 
fic Genesys 20 spectrophotometer (Waltham, MA, USA). 

For determining glucose, sorbitol and acetate, the sam- 
ples were centrifuged for 5 min at 14000 rpm with a 
Sigma 1-15 Microfuge (Osterode am Harz, Germany) and 
the supernatants were used for the analysis. Concentra- 
tions of glucose, sorbitol and acetate in the respective 



supernatants were determined using a Dionex HPLC 
(Dionex, Sunnyvale, USA) with an Organic Acid-Resin 
300 x 8 mm (CSChromatographie, Langerwehe, Germany) 
and a Skodex RI-71 detector. Sulfuric acid in a concentration 
of 5 mM was used as solvent at a flow rate of 0.6 ml/min 
and a temperature of 60°C. 

Software 

OTRs, the calibration factor (K) and simulation examples 
were calculated with MATLAB R2010b (The Math Works, 
MA, USA). The mathematical optimization problems were 
solved using the trust-region-reflective algorithm. Ordinary 
differential equations where solved using the trapezoidal rule. 

Appendix 

Model used for calculating the data of Figure 4: 
Biomass equation: 



dX 
dt 



(A-X 



Substrate equation: 
dS _ 1 

Dissolved oxygen equation: 

dcm 1 

Oxygen partial pressure equation: 



dp { 



'02 



- 02, m 



V - V 



t'Po2 



dt 

Growth rate: 

fi = ft r 



Vg 



OTR 



VlR'T 
Vg 



C02 



(C.I) 



(C.2) 



(C.3) 



(C.4) 



(C5) 



S + K s c 0 2 + K02 
Oxygen transfer rate: 

OTR = k L a-(c 02 * -C02) (C.6) 
Dissolved oxygen at the gas-liquid interface: 

C02* =L 0 2-P02 (C.7) 

With: Oxygen solubility L Q2 = 0.0011 mol/L/bar, mass 
transfer coefficient k L a = 0.082 1/s, volumetric flow into 

the flask V in = 0.67 L/h, volumetric flow out of the flask 

V out = 0.7 L/h, maximal growth rate (i max = 0.26 l/h, half 
velocity constant of the substrate K s = 4 g/L, half velocity 
constant of oxygen K 0 2 = 8 TO" 10 mol/L, oxygen partial 
pressure of inlet po2,in = 0.2095 bar, yield coefficient of the 
substrate Y X s = 0.5 g/g, liquid volume V L = 10 mL, gas 
volume V G = 270 mL, gas constant R = 0.0831 L*bar/mol/K, 
yield coefficient of oxygen Y X Q2 = 41.6 g/mol, temperature 



Hansen et al. Journal of Biological Engineering 2012, 6:1 1 
http://www.jbioleng.0rg/content/6/l/l 1 



Page 11 of 12 



T = 303.15 K, initial conditions: biomass concentration 
X o = 0.5 g/L, substrate concentration S 0 = 45 g/L, oxygen 
concentration in the liquid c 0 2,o = 0.22 mmol/L, oxygen 
partial pressure p 0 2,o = 0.2006 bar. 

Nomenclature 

a, b, c: Polynomial coefficients; C02: Dissolved oxygen in the liquid (mmol/L); 
c 02i o. Initial dissolved oxygen in the liquid (mmol/L)c 0 2: Dissolved oxygen at 
the gas-liquid interface (mmol/L); f: Smoothing splinek; L a: Mass transfer 
coefficient (1/s); K: Calibration factor (bar/mV); K02: half velocity constant for 
oxygen (mol/L); K s : Half velocity constant for the substrate (g/L); L02: Oxygen 
solubility (mol/L/bar); OTR: Oxygen transfer rate (mol/L/h); p: Regularization 
parameter (-); p amb : Ambient pressure (bar); p H20 : Water vapor partial 
pressure in the headspace (bar); P02: Oxygen partial pressure in the 
headspace (bar); p Q 2,o: Initial oxygen partial pressure in the headspace (bar); 
Po2,i n : Oxygen partial pressure of inlet flow (bar); p 0 2,reah Steady state oxygen 
partial pressure (bar); R: Gas constant (bar L/mol/K); RQ: Respiratory quotient 
(-); S: Substrate concentration (g/L); S 0 : Initial substrate concentration (g/L); 
T: Temperature (K); U02: Oxygen sensor signal (mol/L/h); V G : Gas volume (L); 
V L : Liquid volume (L); V m : Molar volume (L/mol); V jn : Volumetric flow into the 
flask (L/h); V out : Volumetric flow out of the flask (L/h); X: Biomass 
concentration (g/L); X 0 : Initial biomass concentration (g/L); Y X0 2: Yield 
coefficient for oxygen (g/mol); Y X s: Yield coefficient for the substrate (g/g) 
uGrowth rate (1/s); u max : Maximal growth rate (1/s); t: Lag time of oxygen 
sensor (h). 
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